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Abstract 

We demonstrate how to match pre-equilibrium dynamics of a 0+1 dimensional quark gluon 
plasma to 2nd-order viscous hydrodynamical evolution. The matching allows us to specify the 
initial values of the energy density and shear tensor at the initial time of hydrodynamical evolu- 
P^H ! tion as a function of the lifetime of the pre-equilibrium period. We compare two models for the 

£Nj | pre-equilibrium quark-gluon plasma, longitudinal free streaming and collisionally-broadened longi- 

tudinal expansion, and present analytic formulas which can be used to fix the necessary components 
of the energy-momentum tensor. The resulting dynamical models can be used to assess the effect 
of pre-equilibrium dynamics on quark-gluon plasma observables. Additionally, we investigate the 
dependence of entropy production on pre-equilibrium dynamics and discuss the limitations of the 
^ c*j standard definitions of the non-equilibrium entropy. 

^ ! PACS numbers: 24.10.Nz, 25.75.-q, 12.38.Mh, 02.30. Jr 
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I. INTRODUCTION 



The goal of relativistic heavy-ion experiments is to produce and characterize the ther- 
modynamical and transport properties of matter at extremely high temperatures. In such 
collisions experimentalists expect to produce a new high-temperature state of matter that is 
a deconfmed plasma of quarks and gluons (QGP). Data from the Relativistic Heavy Ion Col- 
lider (RHIC) are consistent with the formation of a thermalized QGP that exhibits strong 
transverse collective flow. However, there are still many unresolved issues. One of them is 
to determine if the formed plasma is (nearly) isotropic and thermal at early times r ~ 1 — 2 
fm/c. Hydrodynamics predictions for collective flow of the matter are consistent with RHIC 
data using thermalization times in the range r ~ 0.5 — 2 fm/c [l|-[5j]; however, many out- 
standing questions remain. One of the chief uncertainties in hydrodynamical modeling is the 
proper initial conditions to use when integrating the resulting partial differential equations. 
The initial conditions required are the energy density profile, 8 , the components of the fluid 
four-velocity, w M , and the stress tensor Yi' 11 ' at an initial time Thydro- 

In this work we demonstrate how to determine the initial energy density £ and shear n 
in a 0+1-dimensional model. We introduce a pre-equilibrium period in which the system 
develops a local momentum-space anisotropy owing to the longitudinal expansion of the 
matter. After this period we evolve the system using second-order viscous hydrodynamics 
with initial conditions consistent with the pre-equilibrium evolution of the matter. To frame 
the discussion we introduce two proper time scales: (1) the parton formation time, r , which 
is the time after which coherence effects in the nuclear wave function for the hadrons can 
be ignored and partons can be thought of as liberated; and (2) the time at which modeling 
of the system using viscous hydrodynamics, r hydro begins. During the pre-equilibrium stage, 
t < r < T hydro , longitudinal expansion of the matter along the beam axis makes the system 
colder along the longitudinal direction than in the transverse direction, (p 2 L ) < (p\) 0] 
corresponding to a nonvanishing plasma shear n. 

This paper extends previous work in which we introduced pre-equilibrium interpolating 
models [7H9] to calculate the dependence of high-energy dilepton and photon production on 
the plasma isotropization time 0-15]. In the previous analyses the pre-equilibrium stage 



was matched at late times to isotropic ideal hydrodynamical expansion. Here we show how 
to determine the shear n and energy density £ at a proper-time Thydro given a model for the 
evolution of the microscopic anisotropy of the plasma, £ = \(pt) / (Pl) ~ 1 5 where px and 
Pl are the transverse and longitudinal momenta of the particles in the plasma, respectively. 
This is done by matching to the corresponding pressure anisotropy, A = Vt/Vl — 1, and 
energy density, £. Once this matching is performed one can solve the second-order viscous 
hydrodynamical differential equations to determine the further time evolution of the system. 

In Fig. [1] we show the time evolution of the pressure anisotropy and energy density 
assuming Thydro — 1 fm/c resulting from the models described herein. As shown in Fig. [1] 
(left) the magnitude of A is larger in the weakly coupled case starting from the same initial 
pressure anisotropy at r = 1 fm/c. Figured] (right) shows the typical time evolution of 
the energy density using our matching. As shown in this figure in the weakly coupled case 
the 0+1-dimensional plasma lifetime is increased owing to the larger shear viscosity. In the 
body of the text we show how such models are derived and how, specifically, the matching 
at rh y dro is performed. The resulting models can be used as input to predict the effect of the 
pre-equilibrium period on QGP observables such as dilepton and photon production, heavy 
quark screening, etc. 
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FIG. 1: Time evolution of the pressure anisotropy A (left) and energy density (right) for the case 
of collisionally broadened pre-equilibrium dynamics. In this plot we have fixed Thydro = 1 fm/c. 
The dashed line (red) is the case of weakly coupled transport coefficients and the solid (blue) line 
is the case of strongly coupled transport coefficients. 

The work is organized as follows. In Sec. [Til we introduce the second-order viscous 
hydrodynamics differential equations. In Sec lIIH we introduce the two models for the pre- 
equilibrium evolution and discuss how to determine the initial conditions necessary to in- 
tegrate the viscous hydrodynamics differential equations. In Sec. |V] making use of the 
resulting dynamical evolution, we calculate entropy production as a function of r hydro . In 
Sec. IVIj we present our conclusions and outlook. 



II. 0+1 VISCOUS HYDRODYNAMICS 



We consider a 0+1-dimensional system expanding in a boost-invariant manner along the 
longitudinal (beam-line) direction with a uniform energy density in the transverse plane. 
In terms of proper time, r = \Jt 2 — z 2 , and space-time rapidity, ( = arctanh(z/t), the 
second-order viscous hydrodynamic equations are given by [16|, 
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(2.1a) 
(2.1b) 



where £ is the fluid energy density, V is the fluid pressure, EE = IT^ is the (( component 
of the fluid shear tensor, i] is the fluid shear viscosity, t w is the shear relaxation time, and 
Ai is a coefficient that arises in complete second-order viscous hydrodynamical equations in 



either the strongly [13, ll8| or the weakly coupled limit [16|, |19|, |20[ . 



To solve these coupled differential equations it is necessary to specify the initial conditions 
and the equation of state that relates the energy density and the pressure through V = V{£\ 
For 0+1-dimensional viscous hydrodynamics, one must specify the energy density and II at 
the initial time, that is, £h y dro = £( r hydro) and n hydro = n(r hydro ), where r hydro is the proper 
time at which one begins to solve the differential equations. 
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Transport coefficient 


Weakly-coupled QCD 


Strongly-coupled TV = 4 SYM 


fj = rj/S 
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6f]/T 
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(4.1 -> 5.2) rfS/T 


2 rf<S/T 



TABLE I: Typical values of the transport coefficients for a weakly-coupled QGP 19, 3(], 31] and a 
strongly coupled M = 4 SYM plasma [17], [18 ] . 



In the following analysis we assume an ideal equation of state, in which case we have 
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where for quantum chromodynamics with N c colors and Nf quark flavors, Adof = 2(Af? — 
1) + 7N c Nf/2, which, for N c = 3 and Nf = 2 is Adof = 37. The general method used here, 
however, can easily be extended to a more realistic equation of state. 

In the conformal limit the trace of the four- dimensional stress tensor vanishes requiring 
£ = 3V which, using Eq. (12. 2p . allows us to write compactly 



S = (T/7) 4 , with 7 EE 



30 



TT 2 N, 



dof 
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(2.3) 



Likewise we can simplify the expression for the equilibrium entropy density, S, using the 
thermodynamic relation TS = £ + V to obtain S = 4£/3T or, equivalently, 



S 



±p3/4 
37 ' 



(2.4) 



Note that for a system out of equilibrium the full nonequilibrium entropy is modified com- 
pared to (12 .4p . Using kinetic theory it is possible to show that the entropy current receives 
corrections at second order in gradients |2l|]. In the original approach of Israel and Stew- 
art [22[ 23 1, the non-equilibrium entropy is expanded in a series in deviations from equi- 
librium and higher-order corrections are neglected. The Israel- Stewart (IS) ansatz for the 
nonequilibrium entropy is 



S 



is 

noneq 



s 



(2.5) 



where 02 is an a priori unknown function that determines the importance of second-order 
modifications to the entropy current. The IS ansatz satisfies the second law of thermody- 
namics d^S^ > and for massless particles described by a Boltzmann distribution function 
one finds /?2 = T n / (2rj). Recent analyses have shown that, including all relevant structures in 
the gradient expansion, the non-equilibrium entropy contains additional terms not present 



in the simple IS definition of the non-equilibrium entropy |24j-l29l . 

When solving Eqs. ( 12. lap and (12. lbj) it is important to recognize that the transport - 
coefficients depend on the temperature of the plasma and hence on the proper time. We 



summarize in Table [I] the values of the transport coefficients in the strong 17], [l8(] and 
weak [l^, 3(J 31] coupling limits. We point out that in both cases the transport coefficients do 
not satisfy universal relations and therefore, their values can only be taken as estimates. In 
the weakly coupled case the QCD transport coefficients depend on the renormalization scale. 
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In addition, higher-order corrections to some tr ansp ort coefficients from finite-temperature 



33| . At stro ng coup ling;, it has been shown 



38| . We take the preceding 



perturbation theory show poor convergence [3J 
recently that there are corrections for finite t'Hooft coupling [34 - 
estimates in both coupling limits to get a qualitative understanding of what to expect in 
each regime. 

In Table [H the reader should note that in the strong and weak coupling limit the coef- 
ficients and Ai are proportional to r w oc fj/T and Ai oc fj 2 S/T, respectively. This fact 
suggests that we can parametrize both coefficients as 



T 



-2 



Ai = c Xl r] 



where we have introduced the scaled shear viscosity 

fj = Tj/S . 



(2.6a) 
(2.6b) 

(2.7) 



In our analysis we assume that fj is independent of time. 

The dimensionless numbers fj, c n and C\ 1 carry all of the information about the particular 
coupling limit we are considering. Using the ideal gas equation of state [Eqs. 1 12. 3ft and f ]2.4[) 
the parametrization 112. 6p of r n and Ai can be rewritten in terms of the energy density 8 



r - = ^TTI > ( 2 - 8a ) 

X 1 = ^- 2 c Xl f£ 1 /\ (2.8b) 

We use the following values for the transport coefficients in the case of a weakly coupled 
QGP 

10 9 

V = y > Cn = 6 , c Al = - . (2.9) 

For the strong coupling analysis, we use 

V = ^, c^ = 2(2-log2), c Al =2. (2.10) 



A. Pressure anisotropy 

We introduce the dimensionless parameter A, which measures the pressure anisotropy of 
the fluid as follows 

A = ^-l, (2.11) 

r L 

where V T = (T xx + T™) /2 and V L = T zz = -T{? are the effective transverse and longitudinal 
pressures, respectively. If A = 0, the system is locally isotropic. If — 1 < A < 0, the system 
has a local prolate anisotropy in momentum space and if A > 0, the system has a local 
oblate anisotropy in momentum space. 
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In the 0+1-dimerisiorial model of viscous hydrodynamics one can express the effective 
transverse pressure as Vt = V + H/2 and the effective longitudinal pressure as Vl = V — H. 
Using these definitions for Vt and Vl and the ideal equation of state, we can rewrite Eq. 
f l2~TTj) as 

. . 9 / n(r) \ 
A(T) = 2 Ur)-3n(r) )- ^ 



In different limits we have 



9 U(t) 

lim A(t) «-=^-f , (2.13a) 

n<£ v ; 2 £(t) v ; 

lim A(r)w-1, (2.13b) 

n^-2£/3 

lim A(r) « oo. (2.13c) 



At the initial time r = r hydro , A hydro = A(r = r hydro ) is given by 

9 / Ilhydro 
2 \ ^hydro 3 Ilhydro 



Ahydro = I ( c ^Tr, ) , (2-14) 



or, 



ffhydro Q i O A ^hydro ; (2-15) 

6 + Z hydro 

This expression allows one to find the initial condition for the shear tensor component 
ffhydro knowing Ahydro and £h y dro- This relation is the bridge that allows us to match the 
initial condition for the shear tensor Ilhydro from a pre-equilibrium period of the QGP. The 
precise matching is described in Sect. IHH where we derive the connection between A defined 
in Eq. (12. lip and the £ parameter introduced in Ref. (39| . 



III. 0+1-DIMENSIONAL MODEL FOR A PRE-EQUILIBRIUM QGP 

In this section we present two models for 0+1-dimensional nonequilibrium time evolution 
of the QGP: 0+1-dimensional free-streaming and 0+1-dimensional collisionally broadening 
expansion. In each case below we will be required to specify a proper time dependence of 
the hard-momentum scale, phard, and the microscopic anisotropy parameter, £, introduced 
in Ref. [391. Before proceeding, however, it is useful to note some general relations. 



We assume that any anisotropic distribution function can be obtained by taking an 
arbitrary isotropic distribution function fi SO (p) and stretching or squeezing it along one 
direction in momentum space to obtain an anisotropic distribution. This can be achieved 
with the following parametrization 

/aniso(P,£,Phard) = /iso ( vV + £ (P " n) 2 , Phard) , (3.1) 

where Phard is the hard momentum scale, n is the direction of the anisotropy, 1 and — 1 < £ < 
oo is a parameter that reflects the strength and type of anisotropy. In general, Phard is related 



1 Hereafter, we use n = e z , where e z is a unit vector along the beam- line direction. 
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to the average momentum in the partonic distribution function. The microscopic plasma 
anisotropy parameter £ is related to the average longitudinal and transverse momentum of 
the plasma partons via the relation 0^9|, Eo[ 

t-m- 1 - (a2) 

From this expression, one can see that for an oblate plasma (pj) > ^{pD, then £ > 0. In 
an isotropic plasma one has £ = 0, and in this case, phard can be identified with the plasma 
temperature T. 

We now show how to derive a general formula for the time evolution of the microscopic 
plasma anisotropy £ that allows for a nonvanishing anisotropy of the plasma at the formation 
time followed by subsequent dynamical evolution. This is a straightforward extension of the 
treatment presented in Ref. [8| where it was assumed that the plasma was isotropic at the 
formation time. 

In most phenomenological approaches to QGP dynamics it is assumed that the distri- 
bution function at r ~ ro is isotropic, that is, £(r = tq)=0. There is no clear justification 
for this assumption. In fact, in the simplest form of the Color Glass Condensate (CGC) 
model [4l[ the longitudinal momentum would initially be zero. This configuration corre- 
sponds to an extreme anisotropy with £ (or A) being infinite in the initial state. In the 
CGC framework to generate a nonzero longitudinal pressure it is necessary to include the 
next-to-leading-order corrections to gluon production taking into account the effect of rapid- 
ity fluctuations and full three-dimensional gauge field dynamics. There has been progress 
toward the solution of this problem; however, it is still an open question (for recent advances 



in this area see Refs. [42l. |43| and references therein). We also note that, taking into account 



the finite longitudinal width of the nuclei, studies have shown that it may even be possible 



for the initial plasma anisotropy to be prolate at the formation time [44 . 

Here we assume, quite generally, that the microscopic anisotropy parameter, £, at the 
formation time takes on an arbitrary value between -1 and oo given by £o- The initial 
anisotropy £ can be evaluated using Eq. (13. 2p giving 

fr-HSfc- 1 - (3 - 3) 

Writing the longitudinal momentum as 

(Pi) = (pDo + (Spl) , (3.4) 

and using the fact that, in the case of 0+1 dynamics, the average transverse momentum is 
constant 

(Pt) = (Pt)o, (3.5) 
we can rewrite the general expression for £ given in Eq. ( 13. 2 p as 

e=^T-l- (3-6) 



1 + 



(Sp 2 l ) 



Finally, we can parametrize the time dependence of the plasma's average longitudinal mo- 
mentum squared as 
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Comparing with Eq. fl3.4j) . we obtain 



(Spl) (r^ & 



1 (3.7) 



Inserting this into Eq. (I3.6p . we have 

£(£o,T,T ) = (£o + l)(^r) -i (3.8) 

This expression holds for both of the 0+1- dimensional pre-equilibrium scenarios studied 
in this work. In the case of longitudinal free streaming we have 5 = 2, and in the case 
of collisional broadening we have 5 = 2/3. 2 There are other possibilities for the values 
of this exponent associated with the bending caused by growth of the chromoelectric and 
chromomagnetic fields at early times of the collision. We refer the reader to Sec. Ill of Ref. 
Q for an extended discussion of the time dependence of the scaling coefficient 5. 

In a comoving frame, the energy density and pressure components can be determined by 
evaluating the components of the stress-energy tensor, 

T ^ = Swrf f ^^- (3 - 9) 

Using the ansatz, Eq. ( 13. ip . for the anisotropic distribution function and making an appropi- 
ate change of variables, one can show that the local energy density S and the transverse and 
longitudinal pressures Vt and Vl are 

£(?w = t 00 = i f _L + | £ iso(phard) j (3 . 10a) 



Tl(0 £iso(Phard) 
1 (rpxx _|_ rpyy^ 

3 fi + (e-m(o\ 



v T ( P ^,0 = l(T ££ + Tyy) , (s.iob) 



hard J 



2£V £+1 J 

V L ( Phard ,0 = T zz , (3.10c) 

3 f(Z+l)K(0-r 



n so (p 



hard i 



where "P^Phard) and PJfXPhard) are the isotropic transverse and longitudinal pressures and 
£iso(Phard) is the isotropic energy density. 3 The function 1Z(£) is given by 



no = \ 



1 arctan-\/F 
+ 



(3.11) 



2 Owing to the assumption of no dynamics in the transverse plane, collisional broadening can only increase 
the longitudinal momentum in 0+1 dimensions. 

3 We point out that, in general, one cannot identify "P^° (phard ) , T-'f ° (phard ) and £i SO (phard) with their 
equilibrium counterparts, unless one implements the Landau matching conditions. In Appendix we 
show an explicit example where the Landau matching conditions are implemented. 



S 



and in Eq. (I3.10p it is understood that p har d = Phard(£o, r, r , S) and £ = £(£ , t, t , 5). 

Note that for a conformal system the tracelessness of the stress-energy tensor Tjf=0 im- 
plies £ = 2Vt + Vl- This condition is satisfied by Eqs. (I3.10p for any anisotropic distribution 
function (Eq. 13. ip because for an isotropic conformal state V^° L = £i SO /3. 



A. 0+1-dimensional free streaming expansion 

In the free streaming (f.s.) case, the distribution function is a solution of the collisionless 
Boltzmann equation 

p ■ d x f Ls .(p,x) = , (3.12) 

Consider for simplicity that the one-particle distribution function is isotropic at the forma- 
tion time, t = r . 

=f(^±A\ (3.13) 

\ Phard / 

T=TO \ / 

where pt is the transverse momentum, pi is the longitudinal momentum, and Phard is the 
hard momentum scale at r . The hard momentum scale for particles undergoing 0+1- 
dimensional free streaming expansion is constant in time. In comoving coordinates the 
general solution for free streaming expansion in 0+1-dimensional expansion can be written 

as 



fi,Xn-v) | ^Wl + ^sinh 2 (y-C) I • (3.1 I) 

Phaid V T 

where y is the momentum-space rapidity, r is the proper time, and ( is the space-time 
rapidity. 

Using the free streaming distribution function, Eq. f )3.14p . it is possible to show that the 
average longitudinal and transverse momentum-squared values are given by [6, 8] 

I 

^(PT){. s .ocp 2 hard , (3.15a) 
(pi)f.s. oc pLd ^ • (3.15b) 

Inserting these expressions into the general expression for £ given in Eq. (13.21) one obtains 
£/. s .(r) = (1 + £ )(r/r ) 2 - 1. Therefore, 6 = 2 in Eq. (ESJ). 

One can also determine the temporal evolution of the energy density for the 0+1- 
dimensional free streaming case using Eq. (I3.10p . The resulting temporal evolution of the 
relevant variables is jsf 

£ f . s .(r) = (l + £ )(r/r ) 2 -l, (3.16a) 

Phard(T) = Phard , (3.16b) 

£(r) = ft(£ f .,) (^) 4 - (3.16c) 
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B. 0+1-dimensional collisionally broadening expansion 



In the bottom up scenario |6[, it was shown that, even at early times after the nuclear 
impact, elastic collisions between the liberated partons will cause a broadening of the lon- 
gitudinal momentum of the particles compared to the noninteracting, free-streaming case. 
During the first stage of the bottom-up scenario, when 1 <C Q s t <C a 3 / 2 , the initial hard 
gluons have typical momentum of order Q s and occupation number of order l/oc s . Owing 
to the fact that the system is initially expanding at the speed of light in the longitudinal di- 
rection the gluon number density drops like n g ~ Q z s / (a s Q s r). If there were no interactions, 
this expansion would be equivalent to 0+1-dimensional free streaming and the longitudinal 
momentum pi would scale like 1/r. However, when elastic 2 <-> 2 collisions of hard gluons 
are taken into account [6], the ratio between the longitudinal momentum pi and the typical 
transverse momentum of a hard particle px decreases as 



(pi) 



oc (Q s r)- 2 / 3 . (3.17) 



{P 2 T ) 

This implies that for a collisionally-broadened (c.b.) plasma, £ c .b. = (1 + £o)( r /' r o) 2 ^ 3 ~~ 1) 
implying 5 = 2/3 in Eq. (ETBj) . 

The temporal evolution of £, Phard and £ for the case of 0+1 collisionally-broadened 
expansion is jl] 

ec.b.(T) = (l+eo)(r/r ) 2/3 -l, (3.18a) 

Phard(T) = (Phard)o (W 7 ")^ 9 > (3.18b) 

e{r)=H{U) {^f)'- (3.18c) 



C. Relation between A and £ 

In this section, we derive the relation between the pressure anisotropy parameter, A, 
introduced in Sec. Ill A[ and the microscopic anisotropy parameter, £. Combining Eqs. 
( 13.10b)) and (13.10c)) and using = = £i SO /3, we obtain the following expression for A 

= i«-3)+e((l + 0^-l)". (3.19) 

The evolution of A during the pre-equilibrium stage will depend on the kind of model for 
that stage, that is, either free streaming or collisionally-broadened expansion. For small and 
large values of £ 

limA = ^ + 0(£ 2 ), (3.20a) 
5 

hm A = ^ + G(v^). (3.20b) 
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If one uses Eq. (12.13aj) together with Eq. (I3.20a[) . £ can be related with the shear viscous 
tensor during the viscous regime as 4 

£ = f j + o(n 2 ). (3.2i) 

Note that if one uses the Navier-Stokes value of the shear tensor LIns = 4?//(3r) in the last 



relation, the anisotropy parameter can be expressed as [45 



10 77 

£ns = —j + 0(Tl 2 m ) . (3.22) 
D. Matching the initial conditions 

We now match the general evolution of £ from Eq. (13.81) at an intermediate r hydro and use 
this to fix the initial shear tensor that should be used in the viscous evolution. From Eq. 
(13. 8p . the anisotropy parameter takes a nonvanishing value at r = r hydro , 

£ hydro = C(t = r hydro ) = (1 + £ ) fe^j 5 -1. (3.23) 

Once £h y dro is known, the initial pressure anisotropy Ah ydro and initial energy density £h ydro 
can be determined using Eqs. (13.191) and (I3.10al) respectively. It is then straightforward to 
determine n hydro = Il(r = r hydro ) owing to relation (I2.15p . that is, 

n hy dro = II (r = r hydro ) = - y~ — - j Ahydro • (3.24) 

This expression together with £h y dro gives the full set of initial conditions necessary to solve 
the 0+1-dimensional viscous hydrodynamics equations (12. ip . Note that, by construction, the 
initial conditions do not depend on the particular coupling regime we are interested in. This 
is because at leading order the coupling constant cancels out in the case of a collisionally- 
broadened expansion and in the case of free-streaming it is assumed that there is only free 
expansion. As a result n hy d r o and £h y dro depend only on the type of pre-equilibrium scenario 
considered through the exponent S. 



IV. TEMPORAL EVOLUTION INCLUDING PRE-EQUILIBRIUM DYNAMICS 

In Fig. [2J we show the complete temporal evolution of the pressure anisotropy A(r), 
starting from a pre-equilibrium period and matching at 7"h y dro to viscous hydrodynamical 
evolution. In the plot we show three assumed values of 7h y dro- The initial conditions for 
the strong and weak coupling cases are assumed to be the same in both panels. During the 
pre-equilibrium case r < r < Th y dro, A(r) is determined via its relation to £(t) specified in 
Eq. (I3.19p . In Fig. [2] we have shown the case where £ evolves in the collisionally broadened 
scenario, that is, 5 = 2/3. The matching from pre-equilibrium dynamics to viscous evolution 



An alternative derivation of this relation from kinetic theory is presented in Appendix [C] 
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Weak coupling Strong coupling 

1 4n 




x (fm/c) x (fm/c) 

FIG. 2: Temporal evolution of the pressure anisotropy parameter A for three values of Th y d ro £ 
{0.3,1,3} fm/c. We use To = 0.3 fm/c, the initial temperature at the parton formation time 
Phard = Tq = 0.35 GeV, and an initial value for £o=0. Both plots assume the collisionally-broadened 
scenario, and the transport coefficients during the viscous period correspond to the weak coupling 
(right) and strong (left) coupling regimes. 



occurs at Th y d ro , where, owing to the longitudinal expansion of the plasma, a nonvanishing 
value of £ is generated. Using Eqs. (I3.10a[) . (I3.19p . and (I3.24p we use the value of £(Th y dro) to 
determine the initial values of the energy density and shear necessary for integration of the 
viscous hydrodynamical differential equations. From r hydro < r < r fo , A is determined using 
Eq. (I2.12p . It should be understood that during this period of the evolution the energy 
density and shear are the solutions of the viscous hydrodynamical differential equations, 
Eqs. (12. ip . In both the weak and the strong coupling cases the late-time evolution of A 
is given by the Navier-Stokes solution with 11ns = 4t?/(3t). Also, note that if the pre- 
equilibrium evolution results in an anisotropy that is different from A^s (thydro), then the 
system relaxes to the Navier-Stokes solution within a time of the order of r n . 

As shown in Fig. [2] the initial value of A depends on the assumed matching time. As 
rh y dro increases, Ah ydr o and Ilh y dro increase. If the assumed value of Th y d ro is too large then 
one sees an unreasonably fast relaxation to the Navier-Stokes solution. This is true in the 
collisionally broadened scenario depicted in Fig. [2] and also in the free-streaming scenario, 
5 = 2, which we do not explicitly plot. In the free-streaming scenario the longitudinal 
momentum-space anisotropies generated during the pre-equilibrium period are even larger. 

Another issue that arises is that if the initial shear generated by the pre-equilibrium evo- 
lution becomes too large, it can become comparable to the equilibrium pressure V. If this 
is the case, then it is suspect to apply viscous hydrodynamical evolution. However, this is 
not the only possible way to generate unreasonably large shear. Once the hydrodynamical 
evolution begins it is possible to generate large shear during the integration of the hydro- 
dynamical differential equations. This effect is larger in the weakly-coupling the 
values of f] and T n are approximately 10 and 30 times larger than in the strong-coupling 
case, respectively. This is why in Fig. [2] no large values of A are generated in the case of 
strong coupling, whereas the weak-coupling case there are. One other possibility that arises 
is that the initial value of the shear computed from £ will result in the initial condition 
being "critical", meaning that, when the differential equations are integrated, unphysical 
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behaviors such as negative longitudinal pressures are generated |46|, |47(. In our results, we 
check to see the generated initial conditions are critical and indicate wheter this happens in 
the corresponding results tables. 

The evolution shown in Fig. [2] is typical of the time evolution of A in our model. Of 
course, one can vary the assumed value of £ at the formation time and also consider the free- 
streaming case. For the sake of brevity we do not present plots showing these possibilities, 
as the analytic formulas required, Eqs. (13. 8p . ( 13.1 (Jap . ( I3.19p . and (I3.24p . are simple enough 
for readers to implement on their own. These four equations can be used to generate the 
time evolution of the plasma anisotropy for use in phenomenological applications. In Sec. [V] 
we demonstrate how to use the resulting model and calculate entropy generation using it. 



V. ENTROPY PRODUCTION 



In transport theory, the entropy current is defined as [21 

*"(*) = -/ (^^/(^P)0og[/(^P)]-1}- (5-1) 

Contracting the entropy current with the velocity fluid w M , we obtain the entropy density 
S = u^S^ 1 . Nonequilibrium corrections are usually computed by expanding the distribution 
function around equilibrium j2l[. For the anisotropic distribution function, Eq. (13. ip . the 
entropy density can be calculated analytically in the local fluid rest frame using a change of 
variables, giving 

«%hard, C) = ' 

which, unlike typical expressions for the nonequilibirum entropy, is accurate to all orders to 
£. We note, importantly, that our ansatz, Eq. (13. ip does not fall into the class of distribution 
functions describable using the 14 Grad's ansatz, because when expanded around equilib- 
rium, Eq. (13. ip has momentum-dependent coefficients. Therefore, the entropy production 
from our anisotropic distribution will differ from the 14 Grad's method and IS ansatz (See 
Appendixes IB1 and IC1 for a detailed comparison). 

In both the pre-equilibrium and the viscous hydrodynamical periods we use (15. 2p to 
calculate the percentage entropy generation AS/S . We define 

Aff _ T fo <S(r fo ) - r S(r ) _ S { - S ^ 
S t S(t ) S 

where Sf and So are the entropy per unit rapidity evaluated at r = r fo and r = r , respec- 
tively. 

Note that the two models for pre-equilibrium evolution, free-streaming and collisionally- 
broadened expansions, generate no entropy during the pre-equilbrium period. In the case 
of free-streaming it is obvious that there can be no entropy generation. In the case of 
collisionally broadened expansion there is an implicit assumption that there are no inelastic 
processes. Therefore, in both cases there is no entropy generation. This can be checked 
analytically by using Eqs. (13. 16j) and ( I3.18P and computing either the entropy density or 
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FIG. 3: Entropy percentage as a function of Th y dro for a free streaming pre-equilibrium scenario 
in the weak (left) and strong (right) coupling regimes. We use tq = 0.2 fm/c and the initial 
temperature at parton formation time Tq = 0.35 GeV. 

the number density, in which case one finds that both drop like r _1 0, [8|. 5 Of course, 
these models are an idealization and one expects inelastic processes to contribute to entropy 
production during the pre-equilbrium period in a more realistic model; however, this is 
beyond the scope of the current work. 

The entropy produced during the expansion can be used to constrain nonequilibrium 
models of the QGP j48[. The produced entropy depends on the values of the transport 
coefficients and is sensitive to the assumed value of r]/S. Based on the fact that our pre- 
equilibrium models do not generate entropy, one naively expects that if the viscous hydro- 
dynamical period starts later, then less entropy is produced. However, this is only true if 
during the viscous period we have control over the gradient expansion, that is, |IT/7 :? | <C 1. 
Additionally, for fixed phard, the factor of a/1 + £ in Eq. (15. 2p causes the entropy to de- 
crease monotonically as £ — > oo. The competing effects of phard and £ can cause the naive 
expectation described previously to be violated, as we discuss later. 

In all of the subsequent calculations shown below the parton formation time To is chosen 
to be 0.2 fm/c and the initial temperature, T , at that time is taken as 0.35 GeV. We use 
a freeze-out temperature of Tf Q = 0.16 GeV. The values of the transport coefficients are 
summarized in Eqs. (12. 9p and (12. 10p . In what follows, we report the results for entropy 
production when the initial conditions are fixed at the formation time. 



A. Free streaming model 

Fig. |3] shows the entropy percentage as a function of r hydro in the strong (left panel) 
and weak (right panel) coupling limits when free streaming is used as the pre-equilibrium 



5 This result is found only if one uses the exact expression given by Eq. (|5.2j) . If the IS or 14 Grad's 
expression for the nonequilibrium entropy is used, one will find that, even assuming a free-streaming 
plasma, entropy is generated during the pre-equilibrium evolution. This is obviously incorrect so we use 
Eq. (15. 2[) in all cases. 
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FIG. 4: Entropy percentage as a function of Th y dro for a collisionally broadening pre-equilibrium 
scenario in the weak (left) and strong (right) coupling regime. We use tq = 0.2 fm/c and the initial 
temperature at parton formation time Tq = 0.35 GeV. 



scenario. In both case, after a certain time, which we call r c , there is no entropy generation 
in our model. This is because for Th y d TO > t c the system freezes out while still in the pre- 
equilibrium period of evolution. For the free streaming model this time is approximately 
t c = 3.6 fm/c. 

Fig. E] shows that the entropy production depends on the values of the transport coeffi- 
cients. In the strong coupling case, we see that AS/S increases between 0.5 < r hydro < 3 
fm/c goes to zero after r c . The increase in entropy production owes to the rapidly increasing 
value of IIh y dro in the free-streaming case. In the weak coupling case, there is a similar 
behavior but the effect is less pronounced. Again the increase in entropy production can be 
understood if one considers the values of the initial conditions Ilhydro and £h y dro obtained 
from the pre-equilibrium free streaming expansion. For the free-streaming model, for ex- 
ample, the anisotropy parameter £ increases as r 2 . As a consequence, one can obtain large 
initial values for the shear. In fact, care should be taken because the size of n hydro /'Phydro 
can be of 0(1), making the use of a viscous hydrodynamical description after r hydr o suspect. 
Therefore, it is necessary to check the relative size of II and V at the matching time to assess 
the trustworthiness of the hydrodynamical evolution. The constraint on Thydro is stronger if 
one requires instead a more stringent convergence criterium, ^ V/3, during the hydro- 



dynamical evolution |46|. Here we do not apply this stronger condition but mention it to 



make the reader aware of this caveat. 



B. Collisionally broadened model 

In Fig. HI we show the entropy percentage AS/ Sq in the strong (left panel) or weak (right 
panel) coupling cases as a function of Thydro- In the collisionally-broadened case we find 
r c = 2.2 fm/c, which is shorter than the corresponding time in the case of free streaming. 
This is because in the collisionally broadened (c.b.) case, the energy density decreases more 
quickly (£ cb ~ r -11 / 9 ) than in the free streaming case (£f. s . ~ t -1 )- 

Fig. H] shows that, contrary to the case of free streaming, the entropy percentage as Thydro 
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FIG. 5: Entropy percentage as a function of Thydro for a collisionally broadening pre-equilibrium 
scenario and £o=-0.1, i n the strong (right) and weak (left) coupling regime. We use To = 0.2 fm/c 
and the initial temperature at parton formation time Tq = 0.35 GeV. 



is now a monotonically decreasing function of r hydro . This is because the initial values of the 
anisotropics developed during the pre-equilibrium period are smaller for the collisionally- 
broadened case than the free-streaming case [8] . This result is more physical, as one expects 
that as the hydrodynamical expansion expands later (so larger Thydro) then less entropy is 
produced for a given value of the transport coefficients j48|. We note that we also observe 
that the entropy percentage as a function of Thydro drops more quickly in the weak coupling 
regime compared with the strong coupling case. This is a consequence of Eq. ( 15. 2p . which 
shows that for fixed Phard, as £ increases, less entropy is produced. As pointed out in 
Sec. IIVt owing to the different values of the relaxation time r n in both coupling regimes, 
the anisotropy in momentum space is larger in the weak coupling case than in the strong 
coupling case during the viscous period (see comparison in Fig. [2]). 



C. Including initial anisotropies at the formation time 

In the previous subsections we reported results for the case where there was no 
momentum-space anisotropy at parton formation time, that is, £o — 0. Here we relax 
this assumption. In Fig. El we show the result for entropy production for a prolate initial 
distribution with £ = —0.1 in the collisionally-broadened scenario. In both coupling cases, 
AS/ S decreases as T hy dro increases. Because the initial value of the anisotropy is close to 
zero, generally speaking, the behavior of the entropy percentage is similar to the case where 
there is an isotropic initial state (£o = 0). 

For larger values of £o ; the situation becomes more complicated because we do not have 
control of the size of U/V and the system can become "critical". Therefore, for extreme 
initial anisotropies, it is not possible to determine Thydro based on entropy constraints. In 
Table [UJ, we summarize the bounds on Thy dro obtained by varying £ hi the strong and weak 
coupling regimes when one fixes the initial conditions using collisionally broadening expan- 
sion. "Non determined" indicates cases in which the initial anisotropies were so extreme 
as to cause the system to begin generating negative longitudinal pressure. In those cases 
the viscous hydrodynamics is unreliable. The striking conclusion from Table [TTl is that, in 
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AS/ S < 10% 


AS/ S < 20% 


£o 


Weak coupling 


Strong coupling 


£ 


Weak coupling 


Strong coupling 


-0.5 


Thydro > 0.65 fm/c 


Thydro > 0.75 fm/c 


-0.5 


Non determined 


Non determined 





Thydro > 0.45 fm/c 


Thydro > 0.9 fm/c 





Thydro > 0.3 fm/c 


Thydro > 0.35 fm/c 


10 


Non determined 


Thydro > 0.75 fm/c 


10 


Non determined 


Thydro > 0.65 fm/c 



TABLE II: Bounds on Thydro imposed by requiring either a 10% (left) or 20% (right) bound on 
percentage entropy when considering different values of £o an d transport coefficients. We fix the 
initial conditions through 0+1 collisionally broadening expansion, tq = 0.2 fm/c and T= 350 MeV. 



all cases, the lower bound on Thydro owing solely to entropy considerations is larger in the 
weak coupling case than in the strong coupling case. This is caused by the competing effect 
between increasing anisotropy and dropping temperature in Eq. ( 15. 2p and can be seen in 
the fact that the entropy production decreases more rapidly in the weak coupling panels in 
Figs. HI and EJ We point out that these estimates are lower bounds on the minimum Thydro, 
as our models have no entropy generation during the pre-equilibrium period. In a more 
realistic scenario there would also be entropy generation during the pre-equilibrium period, 
which would add to all curves presented in this section. 

In closing we emphasize that the bounds in Table [III do not factor in the constraint that 
the shear should be small compared to the isotropic pressure. As shown in Ref. [8| requiring 
U/V < 1/3 as a convergence criterion for viscous hydrodynamics and assuming an initially 
isotropic plasma (i.e., LThydro = 0) result in a constraint T hy dro > 5.9T _1 in the case of 
a weakly coupled plasma and Thydro > 0.85 Tq" 1 in the case of a strongly coupled plasma. 
Assuming an initial temperature of 350 MeV this gives Th y d ro > 3.3 fm/c in the case of 
a weakly coupled plasma and Thydro > 0.6 fm/c in the case of a strongly coupled plasma. 
Therefore, assuming an initially isotropic plasma, the convergence constraint can be than 
the constraint implied by entropy production. In general one must compare both constraints 
to determine which results in a stronger condition. 

VI. CONCLUSIONS AND OUTLOOK 

In this paper we have presented a model that allows us to match 0+1 pre-equilibrium dy- 
namics and 0+1 second-order conformal viscous hydrodynamics at a specified proper-time 
Thydro- The pre-equilibrium evolution is modeled by either free-streaming or collisionally- 
broadening expansion. We have derived a relation between the microscopic anisotropy 
parameter £ and the pressure anisotropy of the fluid A. This relation allowed us to de- 
termine the initial conditions for the shear II and energy density £ at the assumed matching 
time. The initial values of £ and II depend on the kind of pre-equilibrium model considered 
and, also, on the interval of time over which pre-equilibrium dynamics is assumed to take 
place. The resulting models can be used to assess the impact of pre-equilibrium dynamics 
on a variety of observables such as photon and dilepton production, heavy-quark transport, 
jet-medium induced electromagnetic radiation, etc. 

As a particular application here we have studied entropy generation as a function of 
Thydro- We have derived an exact expression for the nonequilibrium entropy and then used 
this to determine the percentage entropy generation. We have shown that owing to the 
reduction in entropy by a factor of a/1 + £ compared to the isotropic case, it is possible to 
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have more entropy generation in the strongly coupled case. We have summarized our results 
for entropy generation in Table [III by presenting bounds on Th y dro that result from requiring 
the percentage entropy generation to be less than 10% or less than 20% . 

It is possible to extend these studies to higher dimensions and match all components 
of the energy-momentum tensor and fluid four-velocity. To do this it is necessary to spec- 
ify information about the transverse expansion during the pre-equilibrium period and how 
this impacts the anisotropy at early times. One way to approach this problem is to use 
three dimensional (3D) parton cascade models [49-51] or 3D Boltzmann-Vlasov- Yang- Mills 
simulations [52]. Additionally, in our approach there is no entropy generation during the 
pre-equilibrium phase. This is not necessarily true, as inelastic collisions such as 2 — > 3 are 
necessary for chemical equilibration jfj|, 49, 50|| and therefore, their inclusion will produce 
entropy during the nonequilibrium phase of the QGP. Short of this, one can investigate 
simple analytic models such as 3D free-streaming or 3D collisionally-broadened expansion 
and develop analytic models that can be used to determine the necessary initial conditions 
self-consistently. Work along these lines is currently in progress. 
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Appendix A: Landau matching conditions of the anisotropic distribution function 
during the viscous period 

To determine the isotropic equilibrium energy density S eq (T) from a nonequilibrium 
single-particle distribution function, it is necessary to implement the Landau matching con- 
ditions 

£{T) = UpT^Uv , (Ala) 
u^SN" = , (Alb) 
UpST^Uv = , (Ale) 

where is the energy-momentum tensor computed with the equilibrium distribution func- 
tion f eq (x,p), and 5T^ V involves nonequilibrium corrections to the energy-momentum tensor. 
In the 14 Grad's method this constraint is immediately satisfied by construction. In the case 
of the anisotropic Boltzmann distribution, the first constraint, Eq. OA lap , requires 

/ (2^o ( M • 2 ex P [" VP 2 + / Phard] = J ^^(M-p) 2 exp[-p/T] . (A2) 
Performing the integrals on both sides, we find that 
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Phaxd= (ft(0) 1/4 T. (A3) 

Now we expand the anisotropic distribution function to second order in £, and making use 
of the last expression, we have 6 



,/aniso(P,£,T) = expf-vV +fp2/ phard ] (A4) 



cxp _^(^(o) 1/4 v / r+i^ 

« e- p / T (l + e/(i) + e 2 /( 2) ), 
where we use explicitly p z = pcos9. The functions /m and /( 2 ) are given by 

/(i) = ^(1-3 cos 2 9) (A5a) 
/( 2 ) = [5p - 39T + 30(T - p) cos 2 9 + 45(T + p) cos 4 9] . (A5b) 

Replacing the expansion of the anisotropic distribution until C?(£) in the Landau condition, 
Eq. fjATbl 

u ^ N " = ^ S w? e ' p,T f{1) m 



1 -t roa pit 

-^7^ / / ^^(cos^)p 3 e- p / T (l-3cos 2 i 



(2tt 
= 0. 

To 0(C,), the Landau condition, Eq. ( lAlbl) . is satisfied. Note that this condition is not 
expected to hold at all orders because the particle number density is proportional to the 
entropy density. The only way to hold both the energy density and the particle number 
density is by introducing a chemical potential. 

Expanding the anisotropic distribution to 0(t; 2 ) in the Landau condition (I Aid) , we have 

u„5n„ = J -0j- 3 Pe~ p/T (£/(i) + ef(2)) , (A7) 
where fm and fr 2 ) are given by Eqs. (IA50 . Explicitly, we have 

* / Wf pe ~ PlTfw = w?*f I d3pp2e ~ p/T ^ - 3cos2 ^ (A8) 



(2tt) 2 6T 
0, 



j dpp 4 e- p/T J d(cos0)(l -3 cos 2 



(A9) 



6 For practical purposes we expand until second order in £, as we are considering viscous hydrodynamics 
until second order in gradient expansion. 
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+ 30 (T - p) cos 2 9 + 45(p + T) cos 4 0) 

t2 / roc poo \ 

= i =(8/ dpp 5 e- p/T -40T / rfpp 4 e- p/T 

= 0. 

Therefore, up to second order in £, the anisotropic distribution function satisfies the Landau 
condition, Eq. (lAlcj) . 



Appendix B: Entropy from the 14th Grad's Method 

We can evaluate the entropy from the kinetic theory definition using the 14th Grad's 
approximation for the nonequilibrium distribution function 

f(x,p) = f eq (l + 5f), (Bl) 

The dependence of 5f is assumed to be a function of the hydrodynamic degrees of freedom 
£, V, m m , g^ u , and expanded in a Taylor series 

Sfix^jf) = e + + e^p v + 0(p 3 ). (B2) 

By demanding that 5f vanishes in equilibrium, one finds that in the Landau frame and 
assuming massless Boltzmann particles, 5f is given by [22| 

P") = 2T*(£(T)+V(T)) WV + ° ip3) (B3) 
Expanding the expression of the entropy 

S = ~ J 0s / eq (l + 5f) [log(/ cq (l + Sf)) - 1] (B4) 

~ «5(o) + 5(i) + 5(2). 



where 



/Jo 
^ /eq [log(/ eq )-l], (B5a) 



<%) = " / 7^/e q 5/log[/ eq ], (B5b) 



s ^=-\S 0) M5ff - (B5c) 

After replacing the 14 Grad's ansatz in the last expressions, these integrals can be calcu- 
lated analytically if one rewrites them as moments of the equilibrium distribution function 
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feq(x^,p^) = e p / T . After a lengthy calculation, we find: 



1 / 3 T 4 T 4 \ 1 
5(0 = y (— + V2 ) - y (m + V{T)) (B6a) 

= (B6b) 
3 7T 

' 5(2) = lT%j (B6C) 
Using the ideal equation of state, the nonequilibrium entropy is 

4 3 7T v 'K ilV 

Comparing this expression with the IS ansatz for the nonequilibrium entropy, Eq. (12. 5p . we 
find a well-known result for a Boltzmann gas 

For 0+1-dimensional case, where = diag (0, 11/2, 11/2, —II), the nonequilibrium entropy, 
Eq. (JB7]), is 

5„one q = ^(T)-^n 2 . (B9) 

Appendix C: Entropy from the anisotropic distribution ansatz 

Using the kinetic theory framework, we can calculate the entropy from the anisotropic 
distribution ansatz, Eq. (I3.ip . expanding in a Taylor series in terms of the anisotropy 
parameter £ (Eqs. lA~5j) . Since the Landau matching conditions are satisfied by the anisotropic 
distribution function, one can write the shear tensor to first order in £ as 

n " = Sj^k rf! " sf (ci) 

where we use explicitly the first-order correction to the anisotropic distribution function, 
Eq. ( lA5al) . We are interested in the 0+1-dimensional case, where there is just one indepen- 
dent component of the shear tensor II 22 = —II. Calculating the zz component from the last 
expression, we have 



L [ «=» e -p/Tfi _ 3 



tt" = -II = ^ / 7^pie- p/1 (l-3cos^), (C2) 



(2tt) 2 6T 

8 T 4 
~15^' 



P7T roo 

/ d(cos8) cos 2 (1 — 3 cos 2 9) / dpp 



4 e - P /T 
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This expression coincides with Eq. (13.2ip . 



Expanding the entropy to second order in £: 

7 3 



/Jo 
72^/ eq [i + + eh)] [io g [/ eq (i + e/d) + e 2 /( 2) )] - 1] , (C3) 



where 



5(0) + ^(i) + 5(2). 



5(0) = -/ ^/e q [l0g/ eq -l], (C4) 
(^s/w/eqlog/eq, (C5) 



rf3 P , ( (■/"(!) 



\2 



5(2) = y ^/e q (*^- + /(2) log / cq J . (C6) 

After replacing and /(2) (Eqs. ( 1A5I) ) in the last expressions, and with / eq = e _P//T , we 
have 

1 /3T 4 T 4 \ 1 
5 (0 = + ^J^(£(T)+P(T)) (C7) 

5(i) = (C8) 



2 T 



i 



= -T^^ 2 (C9) 



7T 



15 T 

2 ^CO^ 



45 T 

Using the ideal equation of state, the nonequilibrium entropy can be written as 

5 110 n Cq = - j, - — j^"^ 2 (CIO) 



4g(T)/ e\ 

3 t V soy 

4 £(T) f 135 / II \ 2 \ 

3 T V ~ 128 V^M/ / ' 



If one compares this result with the IS ansatz for the nonequilibrium entropy in the tri- 
dimensional case, Eq. (IB9p . the term /3 2 for the anisotropic distribution ansatz can be fixed 
as 

k = §^ (Cll) 

Note that the difference in the values between Eqs. (IB8P and Eq. fIClip comes from the fact 
that the anisotropic distribution is incompatible with the 14 Grad's ansatz. This is because 
at small £ the linear term coming from the ansatz, Eq. (13. ip . cannot be expressed in the form 
a^P/JI^Pis, with a 42 being a momentum-independent quantity. In the case of Eq. (13.11) the 
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corresponding coefficient of p^U^p^ is momentum dependent, and hence Eq. (13. ip does not 
fall into the class of distribution functions describable using the 14 Grad's ansatz. 
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